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The charged particle multiplicities of two- and three-jet events from the reaction 
e+e~ — > Z" — > hadrons are measured for decays to light quark (uds) flavors, 
using the data sample of the OPAL Collaboration at LEP. Using recent theoretical 
expressions to account for biases from event selection, results corresponding to 
unbiased gluon jets are extracted. The unbiased gluon jet data are compared 
to corresponding results for quark jets. We determine the ratio r = Ng/Nq of 
multiplicities between gluon and quark jets as a function of energy scale. We 
also determine the ratio of slopes, r^^^ = (dNg/d3/)/(dNq/dj;), and of curvatures, 
r(2) = (d^Ng/dj/2)/(d2Nq/dj/2), where y specifies the energy scale. At 30 GeV, 
we find r = 1.422 ± 0.051, = 1.761 ± 0.071 and r'^) = 1.98 ± 0.13, where the 
uncertainties are the statistical and systematic terms added in quadrature. These 
results are in general agreement with theoretical predictions. 



1 Introduction 

The mean charged particle muhiphcity of a gluon jet has often been measured 
in the annihilation of an electron and positron to hadrons, e+e^ hadrons. 
The usual method is to select three-jet quark-antiquark-gluon qqg final states 
for which the events and individual jets are defined using a jet algorithm. 
The particle multiplicity of a jet determined with this technique is found 
to depend on which algorithm is employed. Therefore, these jets and the 
associated qqg events are called "biased." In contrast, theoretical calculations 
usually define gluon jet multiplicity inclusively, by the particles in hemispheres 
of gluon-gluon {gg) systems in an overall color singlet. Quark jets are defined 
analogously as hemispheres of quark-antiquark (gg) systems. The hemisphere 
definition of jets yields results which arc independent of a jet finder. Therefore, 
these jets are called "unbiased." Uiibiased gluon jet multiplicity has so far 
been measured only in T Idfl and iP eIqu decays, corresponding to jet energies 
of about 5 and 40 GeV, respectively. 

It is of interest to measure the multiplicity of unbiased gluon jets ai. other 
scales. Such measurements would allow a test of recent predictions W\ from 
Quantum Chromodynamics (QCD) for the scale dependence of the multiplic- 
ity ratio r = Ng/Ng between gluon and quark jets and for the related ratios 
of slopes = (dNg/d?/) and of curvatures r^^^ = (d^Ng/dy^), where Ng and 
Nq are the mean particle multiplicities of gluon and quark jets, y = In [Q / h), 



Q is the jet energy and A is the QCD scale parameter. Recently, a method to 
extract the miiliiplicity of unbiased gluon jets_from biased e^e~ — > qqg events 
was proposed eIIj, extending earlier formalism 113. By combining measurements 
of Ng found from this method with the unbiased measurements for Ng, the 
ratios r, r^^^ and r^^^ can be determined at a variety of scales and used to test 
the corresponding QCD results in a quantitative manner. 

Note that the DELPHI Collaboration has alsfl jjecentl^ presented prelim- 
inary results testing the theoretical formalism of U% (see M) . 



2 Analysis method 

Analytic expressions for the mean particle multiplicity of e'^e^ three-jet 
events, valid to the next-to-leading order ai perturbation theory (NLO, also 
called MLLA), were recently presented in 

^qqg = (L, k^Xu) + \ Ngg (fcj.,L„) , (1) 

Nggg = N,, [Lqq, k^^u) + ^ ^gg ik±,Le) ■ (2) 

The reason for the two different expressions is that there is an ambiguity in 
the definition of the gluon jet transverse momentum when the gluon radiation 
is hard. fc_L and k^ i^^, are the transverse momenta of the gluon with respect 
to the quark-antiquari_system using the definition of cither the Lund (k^^u) □ 
or Leningrad {k±^Lf,) Ii3 groups, ^qqg is the particle multiplicity of a three-jet 
event selected using a jet algorithm. N^^ and Ngg are the multiplicities of 
two-jet qq and gg systems, given by twic<i.Ng and Ng, respectively. The scales 
L, Lqq, k±xu and k±,Le are defined in ffl. 

The multiplicity of the gg system in this formalism, Ngg, depends only 
on a single scale: fc_L,L„ in eq. (jl]) or k±xe in eq. (|^). This dependence on a 
single scale is a statement that Ngg is unbiased, i.e. Ngg {kj_) in eq. or 
is equivalent to the inclusive multiplicity of a gg event from a color singlet 
source produced at the same scale k±, to NLO accuracy. 

The theoretical formalism is based on massless quarks. Therefore, we 
select light quark (u, d and s) events for our study. Three-jet events are defined 
using a jet algorithm. For our standard analysis we employ the Durham jet 
finder. As a systematic check, we use the Cambridge and Lucius jet finders. 
The resolution scale of the Durham jet finder, j/cut, is adjusted separately 
for each tagged uds event so that exactly three jets are reconstructed. For 
the Cambridge jet finder, the resolution scale is again ycut- For Lucius the 
corresponding parameter is djoin- The jets are ordered from 1 to 3 such that 



jet 1 has the highest energy. The angle opposite jet 1 is called 6i, etc. Events 
are retained if the angles between the highest energj^ jet and the other two are 
the same to within 3°, the so-called "Y events" lijilj. For Y events, the three- 
jet event multiplicity and scales Lgg, k± ]^u and k^ ^^, depend only on i?c.m. 
and one inter-jet angle, conveniently chosen to be 9i. For 35° < 6*1 < 120°, the 
range of 9i we employ for our gluon jet analysis, 22 365 events are selected: 
this is our final event sample. For simplicity, we identify the gluon jet by 
assuming it is the lowest energy jet in an event, i.e. jet 3. 

To find the N,, (L, fc_L,L„) terms in eq. (|^), we employ two methods. First, 
for the standard analysis, we perform a direct measurement. Specifically we 
determine the particle multiplicity of two-jet uds events from TP decays as a 
function of the jet resolution scale fc_L,L«- Second, as a systematic check, we 
evaluate the following analytic expression from eI: 

N,^ (J, = N,^ {,]') + (J - J') '^^fjP , (3) 

where J' = k_\_^Lu + Cg with Cq = 3/2. The unbiased term Ng^ {K) is equivalent 
to the mean multiplicity of inclusive e"'"e~ — > hadrons events as a function 
of K — Ec.m.- To find the {Lgg, k±^Lu) terms in eq. M ) , we utilize only 
the second of these methods, i.e. the one based on eq. beacuse a direct 
measurement is not straightforward in this case. 

3 Results 

The leftmost set of plots in Fig. ^ shows our results for the unbiased charged 
particle multiplicities of gg events, N^^-. The solid points in the top plot are 
obtained from eq. (|l|) using the direct measurements of N^^- {L,k±^Lu)- The 
asterisks and open symbols in the bottom plot show the corresponding results 
from eqs. and (||) using the calculated expressions for N^^- (J, k±^Lu) from 
eq. (H). The two sets of results from eq. (solid points in the top left plot 
and asterisks in the bottom left plot) are seen to be very similar. 

Included in the leftmost plots_are direct measurements of mibiascd gluon 
jet multiplicity from the CLEO uu and OPAL Collaborations eI. Also shown 
are the predictions of the Herwig and Jetset Monte Carlo event generators 
for the inclusive charged particle multiplicity of gg events. It is seen that 
the Monte Carlo predictions describe the direct measurements of unbiased 
gluon jet multiplicity (triangle symbols) well. The results from the present 
analysis based on eq. (|l]) are also well described by the Monte Carlo curves. In 
contrast, the results from eq. (§) (open symbols in the bottom leftmost plot) 
are generally well above the Monte Carlo predictions, and - if extrapolated to 
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Figure 1. Leftmost two plots: the mean charged particle multiplicity of unbiased gg events 
as a function of scale. Central plot: the ratio of slopes r between unbiased gluon and uds 
quark jets. Rightmost two plots: the corresponding results for the ratios of slopes r^^^ and 
of curvatures r'^). 



lower and higher energies - appear inconsistent with the direct measurements 
from CLEO and OPAL as weU. We therefore conclude that the equation based 
on the Lund definition of the gluon jet scale, eq. (|l|), yields results which are 
more consistent with other studies than the equation based on the Leningrad 
definition, eq. (||). We henceforth restrict our analysis of gluon jets to the 
former set of results. 

We employ the following procedure to determine r, and r^^^ from 
experiment. The corrected unbiased quark and gluon jet multiplicities are 
separately fitted using thc.pext-to-next-to-next-to-leading order (3NL0) ex- 
pressions for Ng and Ng fl^til, respectively. To better constrain the results, 
the direct measurements of unbiased gluon jet multiplicity from CLEO and 
OPAL are included in the gluon jet fit. The ratio of the fitted expressions 
for Ng and defines r. We calculate the first and second derivitives of the 
analytic equations for Ng and Ng with respect to y and evaluate the resulting 
expressions using the corresponding fitted parameter values. The ratios of 
these terms define r*^^^ and r^^\ 

Our results for r are shown in the central plot in Fig. |l|, those for r^^^ 
and r*^^) in the two rightmost plots. The central results are indicated by solid 
curves. The shaded bands show the statistical uncertainties. The overall 
uncertainties, with statistical and systematic terms added in quadrature, are 
shown by the open bands. At 30 GeV, a typical scale in our analysis, we find 
r = 1.422±0.006±0.051, r^^) = 1.761±0.013±D.070 and r^^) = 1.98±0.02±0.13. 
This is consistent with the QCD prediction Q that r < r^^) < r'^' < Ca/Cf = 



2.25 for the scales accessible in our study, |-| 

The 3NL0 predictions of Capella et al. 13 for r, and are shown 
by the long-dashed curves in the central and rightmost plots of Fig. |l|. The 
analytic predictions for r and r^^^ exceed the corresponding experimental re- 
sults by about 22% and 6%, while the theory agrees with the data for 
A second QCD prediction for r versus scale was recently presented by Lu- 
pia and Ochs The predictions of this calculation for r and r*^^^ are shown 
by the short-dashed-dotted curves. At 30 GeV, r is predicted to be 1.45, in 
agreement with our measurement at that scale. The corresponding result for 

is 1.64, about 7% below the data. Finally, theoretical predictions for x, 
r^^-* and A"^^ can be derived from the formalism of Eden and Gustafson u. 
These predictions are shown by the long-dash-dotted curves. The results are 
seen to be in good overall agreement with the data. We note, however, that 
these predictions are based on the experimental measurements of quark jet 
multiplicities. Therefore, the predictions we derive based onO are not entirely 
independent of the data. 

In conclusion, we find overall agreement between the experimental and 
theoretical results. 
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